Fluctuations of temperature gradients in turbulent thermal convection 
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Broad theoretical arguments are proposed to show, formally, that the magnitude G of the temper- 
ature gradients in turbulent thermal convection at high Rayleigh numbers obeys the same advection- 
diffusion equation that governs the temperature fluctuation T, except that the velocity field in the 
new equation is substantially smoothed. This smoothed field leads to a —1 scaling of the spectrum 
of G in the same range of scales for which the spectral exponent of T lies between —7/5 and —5/3. 
This result is confirmed by measurements in a confined container with cryogenic helium gas as the 
working fluid for Rayleigh number Ra = 1.5 x 10 . Also confirmed is the logarithmic form of the 
autocorrelation function of G. The anomalous scaling of dissipation-like quantities of T and G are 
identical in the inertial range, showing that the analogy between the two fields is quite deep. 

PACS numbers: 47.27.Te; 47.27.Jv 
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While statistical properties of temperature fluctuations 
in turbulent Rayleigh-Benard convection have received 
considerable experimental and theoretical consideration 
(see, for instance, and the references therein), 

corresponding properties of temperature gradients are 
still unexplored from both theoretical and experimental 
points of view. In the present paper, we study statisti- 
cal properties of temperature gradients emphasizing their 
qualitative and quantitative similarity to those of the 
temperature fluctuations themselves. Theoretical consid- 
erations will be based on an equation to be derived for 
the magnitude of temperature gradients, and the results 
will be compared with measurements in turbulent con- 
vection in a confined container of circular cross-section. 
We would like to emphasize crucial difference between 
spectrum of the temperature gradient and spectrum of 
absolute value (magnitude) of the temperature gradient. 
The former one can be readily estimated using the Tay- 
lor's hypothesis. Namely, the Fourier transform of the 
time series of the temperature derivative is simply pro- 
portional to omega multiplied by the transform of the 
temperature time series. If this is the case, there should 
simply be a factor omega squared between the two spec- 
tra. If, however, we consider spectrum of absolute value 
of the temperature gradient, then there is no straight- 
forward relation between this spectrum and spectrum of 
the original temperature time series. In the last case one 
need use of physics, and the spectrum (autocorrelation 
function) of the absolute value of the temperature gra- 
dient can give additional information about the thermal 
convection process. 

The measurement apparatus has unity aspect ratio. 
The sidewalls of the apparatus are insulated, and the 
bottom and top walls are maintained at constant tem- 
peratures; the bottom wall is held at a slightly higher 
temperature A than the top wall. The working fluid 
is cryogenic helium gas. We measure temperature fluc- 
tuations at various Rayleigh numbers towards the up- 
per end of this range, in which the convective motion 



is turbulent, but use the data obtained at the Rayleigh 
number Ra = 1.5 x 10 11 in the present paper. Time 
traces of fluctuations are obtained at a distance of 4.4 cm 
from the sidewall on the center plane of the apparatus. 
This position is outside of the boundary layer region for 
the Rayleigh number considered here. At this Rayleigh 
number, the mean wind (which is the large-scale circula- 
tion within the convection apparatus) is well developed 
so Taylor's hypothesis can be employed when necessary. 
More details of the experimental conditions and measure- 
ment procedure can be found in Ref. 0. 

In thermal convection, the temperature field T(r, t) is 
convected by the velocity field v(r, t), which itself is gen- 
erated by density differences set up between the top and 
bottom walls. We will consider incompressible flow obey- 
ing V • v = (with unit density for simplicity). The 
relevant equations under the Boussinesq approximation 
are 

= -Vp - (v • V)v + t/V 2 v + agTz, (1) 



Of 



dT 

— = _(v V)T + £>V 2 T. 
at 



(2) 



Here p, v, D, a, g and z are, respectively, the pressure, 
kinematic viscosity, thermal diffusivity, isobaric thermal 
expansion coefficient, acceleration due to gravity, and the 
unit vector in the upward direction. Equation (2) is the 
standard scalar advection-diffusion equation, except that 
the velocity field is coupled to the temperature field. This 
'active' nature of the temperature fluctuations in convec- 
tion makes their statistical properties different from those 
of a passive scalar advected by a turbulent velocity with 
no back reaction. Restricting attention, for simplicity, to 
Prandtl numbers of the order unity, the experimentally 
measured spectral density of temperature fluctuations in 
the inertial (convective) range rolls off at a rate that is 
closer to —1.4 than to —5/3, the latter being the case for 
passive scalars in three-dimensional homogeneous turbu- 
lence HiSi0. 



2 



The equation for temperature gradients G = VT can 
be readily derived from (2) as 

^i = - v .^l-^i G . +D ^l ( 3 ) 
dt 3 dxj dxi 3 dx 2 - ' 

with the indices i and j representing the space coordi- 
nates, and the summation over repeated indexes is as- 
sumed. The magnitude G of the temperature gradient is 
determined by G = Gn, where n is the unit vector with 
its direction along vector G. Multiplying both sides of 
Eq. (3) by m, making summation over i, and taking into 
account of the fact that nf = 1, we obtain 

-^j- = — (v • V)G + D\7 2 G - AG, (4) 

which is formally similar to Eq. (2) except for the last 
term in (4). The coefficient A in this term has the form 

dvi fdnA 

Let us now search for circumstances under which the 
last term in Eq. (4) is small in the inertial range. The 
second term in A is assured to be small because the diffu- 
sivity D is small. But the nature of the "stretching" part 
on the right hand side of Eq. (5) is not apparent without 
further considerations. 

As a further step, let us make the following conditional 
average of Eq. (4). Fix the magnitude G in the vector 
field G = Gn while performing the average over all real- 
izations of the direction vector field n permitted by equa- 
tion (3). Let us denote this ensemble average as (...) n . 
From the definition, this averaging procedure does not 
affect G itself, but modifies the velocity field v, which in 
turn modifies the coefficient A in Eq. (4). We may write 

^ = -((v) n -V)G + ^V 2 G+(A) n G. (6) 
at 

It is worth noting that the solutions of Eq. (3) satisfy Eqs. 
(4) and (6), but not all possible formal solutions of the 
Eqs. (4) and (6) satisfy Eq. (3); similarly, not all formal 
solutions of Eq. (6) satisfy Eq. (4) while all solutions of 
Eq. (4) do satisfy Eq. (6). In particular, the solutions 
of Eqs. (4) and (6) are the same only if: (a) the initial 
conditions for the two equations are the same, and (b) if 
realizations of (v}„ and of (A) n , related to these initial 
conditions by the conditional averaging procedure, are 
taken from the applicable solutions of Eq. (4). 

It is difficult to guess a priori when (A) n is negligible, 
because there is no small parameter for the stretching 
part of A. Therefore, let us consider a generic set of con- 
ditions, presumably for the inertial (convective) range, 
which can result in (riinjdvi/ dXj) n — 0. This can be a 
combination of isotropy, which yields 

(n i n : j) n = (i^j) (7) 



and 

<n?)n = (n 2 2 ) n = (nj) n , (8) 
and directional randomness determined by equation 

(rknj<p)n = (ninj) n (ip) n , (9) 

where ip = dvk/dxi for arbitrary k and I. One should not 
mix the conditions (7)- (9) with global isotropy and statis- 
tical independence of the velocity and temperature gradi- 
ents. In particular, (7)-(9) can be satisfied in the inertial 
range even in the presence of strong global anisotropy 
and correlation between gradients (see below for more 
comments) . 

If we use conditions (7)- (9) in the presence of the in- 
compressibility condition dvi/dxi = we obtain 

<*>■— D <(^)'>- (10) 

That is, the formal difference between Eq. (2) for T and 
the conditionally averaged Eq. (6) for G is reduced to the 
"A" term with the A given by Eq. (10). Equation (6) can 
then be reduced, in Lagrangian variables, to 

dF = (A)nG ' 

with the "multiplicative noise" (A) n given by Eq. (10). 

Weak diffusion of Lagrangian "particles" can be de- 
scribed as their wandering around the deterministic tra- 
jectories. The introduction of a weak diffusion is equiv- 
alent to the introduction of additional averaging in Eq. 
(11) over random trajectories [l3T |. The small parameter 
D in (10) and (11) will then determine a slow time in 
comparison with the time scales in the inertial range and 
will therefore not affect scaling properties of G in that 
range. 

We should emphasize that the conditional average in- 
dicated by (. . .)„ is quite different from the global av- 
erage indicated by (...). Because of this, the quantity 
G in Eq. (6) remains a fluctuating variable. To elimi- 
nate the stretching part from the conditionally averaged 
coefficient (A) n , one does not need to satisfy conditions 
(7)-(9) for all realizations of the temperature gradient 
field G, but only for the subset of realizations that gives 
the main statistical contribution to the spectrum of the 
magnitude G in the inertial range. Therefore, conditions 
(7)-(9) could well be violated globally without affecting 
the main conclusion. 

The essential point here is that the conditionally aver- 
aged velocity (v) n is smoothed substantially in compari- 
son with v, while the fluctuation of G itself is still rapid in 
the diffusion-advection equation (6) (because it remains 
in tact under the conditional average, by virtue of its 
definition). Under these typical circumstances, the natu- 
ral expectation (see, for instance, Ref. [14| and references 
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FIG. 1: Spectrum of a one-dimensional surrogate of the mag- 
nitude of the temperature fluctuations gradient in thermal 
convection. The straight line is drawn to indicate the power- 
law spectrum (12). 
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FIG. 2: Autocorrelation function C(r) plotted against logr, 
computed for the same data used for the spectral calculations. 
The vertical arrow here and in Fig. 3 indicates the end of the 
inertial range. 



therein) is that the spectral law for a quantity governed 
by the diffusion-advection equation has a "—1" region. 
The result owes itself to the pioneering work of Batch- 
elor 0], who applied this general idea to the viscous- 
convection range of passive scalar fluctuations. While 
the two contexts are quite different, they are the same in 
the sense that the velocity field is smooth. 

There is another way of deducing the —1 power law. 
We recall that the spectral density fluctuations of tem- 
perature in this region of scales has an approximately 
— 7/5-ths slope. This slope can be derived by the di- 
mensional considerations used by Bolgiano . We may 
apply the same reasoning to Eq. (6). Thus, introduc- 
ing an analogy of the dissipation rate for G, namely 
X* = dG 2 /dt, Bolgiano's dimensional arguments yield 
the scaling law for G to be 



E G (k) ~ <x*> 4/5 (a.<?r 2/5 fc- 



(12) 



in the inertial range. Figure 1 shows the corresponding 
spectrum observed for a one-dimensional derivative of G 
in our experiment. For the effectively smoothed velocity 
field, the space autocorrelation function can be charac- 
terized by a logarithmic behavior given by 



C(r) = (C(r)G(O)) 



r 



(13) 



or, using Taylor's hypothesis, in terms of r by 

C(r) = (G(r)G(0))~ln(^). (14) 

This is seen from Fig. 2 to apply quite precisely for the 
data. In our approach, the active character of the tem- 
perature in the convection manifests itself through non- 
trivial properties of the locally averaged velocity field in 



eq. (6). As a consequence, for instance, the very sig- 
nificant input scale L or tq for the autocorrelation func- 
tion has approximately the same size as typical size of 
the largest plums observed in the convection. The input 
scale can be readily calculated from Fig. 2 (to — lOsec). 

The similarity of the T and G fields can also be seen for 
the scaling of the dissipation rate itself. The local tem- 
perature dissipation can be characterized by a gradient 
measure as 



X{r) 



(15) 



where fl r is a subvolume with space-scale r. The scaling 
law for the moments of this measure, 



(X(r) p ) 
( X (r))P 



(16) 



is an important characteristic of the dissipation intermit- 
tency 0,0]. Using Taylor's hypothesis, one can define 
the local dissipation rate as 



(17) 



and the corresponding scaling of the moments of the local 
dissipation rate as 



(X(r) p ) 
(X(r)>* 



-—Hp 



(18) 



Analogous considerations can be brought to bear for 
the magnitude of the temperature gradient 



X*( r ) 



(19) 



4 



A 




10 100 



under which the equation is identical to that governing 
the temperature itself. The main difference is that the 
velocity appearing the new equation, being a conditional 
average, is a smoothed field. For the advection-diffusion 
equation governed by a smooth velocity field, it is natu- 
ral to expect a power-law spectrum with a slope of — 1; 
measurements of the magnitude of the spectral density 
of G are consistent with this expectation. The correla- 
tion function of G shows a logarithmic behavior, also as 
expected. Finally, the scaling of the square of the deriva- 
tive of G has scaling exponents that are identical to those 
of the temperature itself, confirming that a deep analogy 
exists between T and G in the inertial range. 

We thank J. Schumacher and V. Yakhot for brief dis- 
cussions at an early stage of the work. 



FIG. 3: Normalized moments (x( r ) P ) an d {x*( T ) P } against r 
for the data obtained in convection (p = 2, 3, 4). The straight 
lines drawn to indicate scaling are best fits to the data on the 
left of the vertical arrow. 
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FIG. 4: Intermittency exponents fj, p and fi* extracted as 
slopes of the straight lines in Fig. 3 (equation (18)). 



with exponents fx*. We show in Fig. 3 the dependence 
of the normalized moments (x( T ) p ) an d (x*( T ) p ) 011 T 
calculated for the data obtained in thermal convection. 
The slopes of these straight lines provide us with the val- 
ues of the intermittency exponents \i v and /it*, which are 
shown in Fig. 4. The two sets of intermittency exponents 
obtained for T and G are very close. 

In summary, we have derived an equation for the mag- 
nitude of the temperature gradient G in thermal con- 
vection, and shown that there are general circumstances 
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